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Abstract 


We  consldui*  the  equilibrium  statistical  mechanics  of  classical 
fluids  in  v/hicli  the  potential  energy  Is  decomposable  Into  repulsive 
pair  interactions.  A  generalized  cluster  expansion  Is  derived  relating 
the  thermodynamic  and  structural  properties  of  sucii  systems  to  those 
of  the  hard  spiierc  fluid.  The  expansion  is  ordered  by  a  softness 
parameter  ^  v/hich  is  essentially  the  range  of  intermolecular  distances 
in  v/hich  the  difrcrence  betv/cen  the  Mayer  f  functions  for  the  repulsive 
potential  and  an  appropriate  reference  hard  sphere  potential  is  non¬ 
zero.  The  first  ( lo./est  order)  approximation  generated  by  the  expan¬ 
sion  equates  the  free  energy  and  y( r)  for  the  fluid  to  the  respective 
functions  appropriate  to  a  system  of  hard  spheres  v^ith  diameter  d. 

t 

Here  y(  r)  =  g(  r)  exp[  ■♦3u( r)l ,  v/here  g(r)  and  u(  r)  denote  the  radial 
distribution  function  and  repulsive  pair  potential  respectively.  A 
prescription  is  given  for  choosing  a  temperature  and  density  dependent 
dicmalcsr  d  of  the  reference  hard  sphere  fluid  so  that  the  first  approxi¬ 
mation  for  the  free  energy  contains  errors  of  order  5*  ®nlyt  ond  the 
corrections  to  the  first  approximation  for  g(r)  are  of  order 
The  methol  is  used  to  calculate  the  properties  of  a  fluid  v/hose’ Inter¬ 
molecular  potential  varies  as  r*'*.  The  repulsive  potential  that  produces 
the  repulsive  forces  in  the  Lennord-Jones  potential  is  also  studied. 

Since  the  properties  of  the  hard  sphere  fluid  are  known  from  the  re¬ 
sults  of  computer  calculations  and  conveniently  summarized  by  analytic 
equations,  llic  application  of  the  first  approximation  is  numerically 
very  simple.  With  this  approximation,  the  results  obtained  for  both 
model  systems  agree  closely  with  those  obtained  by  Monte  Carlo  cal¬ 
culations. 


I.  Introduction 


At  least  since  the  time  of  van  der  Waais  and  Maxwell,  the  hard 
sphere  or  "billiard  bail"  model  has  been  thought  to  represent  many 
of  the  essciilial  features  of  realistic  intcrmoiccular  repulsive  forces 
in  equilibrium  and  non-equilibrium  phenomena.  As  a  result,  much  effort 
has  been  devoted  to  understanding  the  properties  of  hard  sphere  fluids. 

The  Percus-Ycv ick  equation*  for  hard  spheres  has  been  solved  exactly,* 
giving  convenient  and  accurate  i nf orn'.a t i on  about  the  equation  of  state 
and  pai r  distribution  function  at  low  densities.  The  scaled  particle 
theory*  is  of  comparable  accuracy  for  the  thermodynamic  prc>perties. 

To  extend  our  knowledge  to  higher  densities,  Monte  Carlo  and  molecular 
dynamics  computer  calculations  have  been  perfonne'd. Recently,  analytic 
expressions  for  the  equilibrium  properties  have  been  presented  which 
accurately  sur.imarize  the  results  of  the  computer  calculations 

at  moderate  and  high  densities  and  v/hich  have  the  proper  functional 
format  low  density.^  In  addition,  molecular  dynamics  calculations  have 
shown  that  hard  spheres  exhibit  a  first  order  phase  transition,* 
which  is  similar  in  many  respects  to  the  liquid-solid  transition  of 
real  fluids,**’ 

The  hard  sphere  model  is  obviously  an  idcol  ization,  and  the  purpose 
of  this  paper  is  to  discuss,  within  the  context  of  equilibrium  statistical 
mechanics,  the  relationship  between  this  idealized  model  and  the 
smoothly  varying  repulsive  forces  found  in  real  fluids.  The  usual 
method  fordoing  this  is  a  perturbation  theory  of  some  type,  and  several 
authors  have  developed  theories  to  relate  hard  sphere  data  to  the  properties 
of  fluids  with  other  repulsive  potentials.  Rcx-/linson  considered  fluids  whose 
repulsive  intermolecular  potential  varied  as  the  inverse  n-th  power  of 


Hu  expnivJod  the  thermodynamic  proper- 


thc  Inturmoluculor  ricparatlon. '  ^ 
ties  In  powers  of  l/n,  the  lowest  order  result  corresponding  to  n»»  which 
Is  equivalent  to  the  hard  sphere  potential.  Barker  and  Henderson* have 
shown  how  to  9cnerall2e  the  Ro’wllnson  method  to  apply  to  a  wide  class 
of  repulsive  potentials  ti^(r).  In  this  generalized  Rcxvllnson  method 
a  tcmperaturc-dopcndcnt  hard  sphere  dianteter  Is  calculated  according 
to 
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and  the  pressure  and  free  energy  arc  set  equal  to  those  of  the  associated 
hard  sphere  fluid  v/lth  the  same  temperature  and  density  and  with  a 
diameter  d  .  Recently,  Darker  and  Henderson  have  modified  this  basic 
approach  by  dividing  the  repulsive  potential  Into  two  parts  at  some 
radius  r^ii.  The  contribution  to  the  thermodynamic  properties  from 
the  portion  r<|i  Is  calculated  by  the  generalized  Rowllnson  method 
while  that  from  the  portion  r>^  Is  given  by  the  first  term  of  the 
high  temperature  expansion.  Finally,  the  parameter  ^  Is  chosen  to 
minimize  the  excess  free  energy  per  particle  as  given  by  this  combination 
of  methods. 

In  our  recent  discussion  of  the  equillbriurr  structure  of  dense 
simple  liquids'*,  we  Introduced  a  new  method  to  relate  both  the  thermo¬ 
dynamic  properties  and  the  equilibrium  structure  of  a  general  repulsive 
system  to  those  of  a  hard  sphere  system.  We  presented  physical  and 
heuristic  arguments  to  Justify  the  procedure.  In  this  paper  we  show 
that  the  approximation  method  Is  In  fact  the  first  step  In  a  rigorous 
and  systematic  procedure.  In  Sec.  II  we  Introduce  and  analyze  the 


theory,  which  Is  boscd  on  o  gnnerni  Izotlon  of  tho  fomMlar  Mayer  cluster 
theory.  In  Sec.  Ill  wo  demonstrate  the  accuracy  of  our  method  for  the 
thermodynamic  properties  of  two  systems  composed  of  "soft  sphere" 
particles.  The  first  system  Is  a  fluid  In  v/hich  the  Intermolecular 
potential  Is  Inversely  proportional  to  the  tv>;elfth  pa-zer  of  the  Inter¬ 
molecular  soporotlon  (the  r"'*  potential).  In  the  second  system,  the 
potential  Is  chosen  to  have  tho  same  repulsive  forces  as  those  of  the 
Leonard- Jones  potential.  In  both  cases,  the  results  obtained  compare 
favorably  with  those  found  by  Monte  Carlo  calculations.  The  paper  Is 
concluded  In  Sec.  IV  with  a  discussion  of  our  method  for  relatInQ  hard 
sphere  and  soft  sphere  fluids. 
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1 1 . _ I  It 

A  convent Iciiol  muIoI  for  •  liquid  Is  o  ilatslcal  system  of  N  Inter- 

octlitij  p.nrticics  In  <i  volunNt  V,  for  which  the  total  poientlol  energy,  U, 

Is  assiirv.'d  to  be  «i  num  of  pair  energies  which  depend  only  on  the  scaler 

disitinccs  between  particles 
N 

U  «  I  u(r|j)  . 

KJ 

I 

When  usln^  the  canonical  ensemble,  the  central  problem  for  this  model  Is 
to  determine  ho.;  the  excess  Helmholtz  free  energy  depends  on  the  potential 
u.  Wc  let 

-  -{lAA/V  -  V“  In  a 
Q-v““/d/  n  v(r,.)  , 

Kj  'J 

1 

p  «  N/V  Is  the  thermodynamic  particle  density,  9*'  I*  BoUzmonn's  con¬ 
stant,  k,  times  the  temperature,  T,  endtp(r)  »  exp(-9u(r)]  Is  the  usuel 

u 

Boltzmann  factor  for  the  pair  potential  u.  Tlie  notation^  denotes  the 
positions  of  the  N  particles,  and  AA  Is  the  excess  Helmholtz  free  energy 
(compared  with  an  Ideal  gas  at  the  same  temperature,  density  and  volume). 
From  ^(p,9;qp)  wo  can  calculate  the  pressure,  entropy,  and  other  therms 
dynamic  properties  by  straightforward  differentiations  with  respect  to 
p  and  9.  Also,  the  radial  distribution  function,  g(r),  can  bo  deter¬ 
mined  by  functional  differentiation  with  respect  tO(p(r): 


(pV2)  g(r)  =q)(r)  (6(?(p,9;q))/B9{»:))  . 
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For  iiiQ  prttsont  purposot,  It  It  foore  convenient  to  deal  not  with  g(r)  directly 
but  rathor  with  y(r),  ««hich  is  defined  at 

y(r;p,a::p) 

Thus,  we  have 

y(r)  ^  (?/p*)  (6^/Mi))  .  (2.3) 

Note  that  y(r)  is  a  more  sla^tly  varying  function  of  r  tlian  is  g.  (In  fact, 
even  for  hard  sptoro  potentials,  y(r)  is  finite  and  continuous  for  ell  r,  as  can  bo 
seen  from  an  inspuction  of  its  virial  series  or  from  the  fact  that  it  is 
aquivalent  to  the  poir  distribution  function  of  cavities  in  a  hard  sphere 
fluid.  •>) 

For  soft  sphere  particles  %«hich  have  continuous,  positive  and  puraiy 
repulsive  potentials,  the  function  (p^(r)  Is  Indicated  scheiMtically  in  Figure 
la.  (in  the  fol lazing  discussion,  the  subscript  s  always  denotes  a  quantity 
for  a  soft  sphere  fluid.)  The  function  rises  smoothly  from  aero  at  smell 
distances  (where  wo  assuao  u^  is  extremely  large)  to  unity  at  large  distances. 

For  comparison,  the  9^  function  for  hard  spheres  viith  diameter  d  is  shown  in 
Flguro  lb.  (Ue  denote  all  quantities  for  a  hard  sphere  fluid  with  a  sub¬ 
script  d.)  If  the  continuous  potential  of  interest  is  harshly  repulsive,  the 
9^  function  is  very  similar  to  the  9^  function,  and  since  ^depends  only  on 
9,  m  might  expect  that  the  thermodynamic  properties  and  y(r)  for 
a  harshly  repulsive  potential  might  bo  quite  similar  to  the  corresponding 
quantities  for  hard  spheres,  provided  that  the  dlamoter  d  is  chosen  properly. 

For  a  reasonable  choice  of  d,  9|(r)  *9^(r)  is  effectively  non-zero  only  for 
a  small  range  of  values  of  r  (See  Figure  Ic. ).  If  that  range  is  $d,  then  ( 
is  «i  dlmonsionless  peramstar  which  Is  a  measure  of  the  softness  of  the  potential. 


Thn  l.iriicr  r  Is,  the  toftur  Is  tlw  potent IaI,  and  6n1y  for  hord  tpharet  Is 
cqudi  to  zora  Wo  will  now  develop  on  expression  for^  In  powers  of  this 
softness  porr.mctcr.  To  do  this  we  need  to  know  hot*/ ^ changes  when  9  Is  changed. 

In  (jenorol,  If  4>(r)  Is  changed  from  one  function  to  another  (I.e.  If 
tp(r)  Is  changed  fron  to  ps(r)  -•  £^r)),  tite  resulting  change  In  Ci  can 

be  expressed  in  terms  of  0  functional  Toylor  series: 

fl{o.P;vo  '  M 


♦  ...  (2.4) 

Those  functional  derivatives  con  be  evaluated  In  several  ways  (see  Appendix 
A).  The  first  derivative  Is  given  In  Eq.  (2.3).  The  second  derivative  Is 

4  (pV<)  ''*'J'*lt**<X'**X*  *<£»«*£’>  y.(r,»)  J,U)(r*) 

(2.5) 

where  6(^  Is  the  three  dimensional  Dirac  delta  function,  and  Je^  ^  and 
are  certain  correlation  functions  (see  Appendix  A).  The  subscript  zero 
Indicates  that  the  correlation  functions  are  those  for  the  fluid  In  which 
(The  p  and  g  dependence  of  the  correlation  functions  In  Eq.  (2.5) 
shoyld  be  understood  even  though  It  Is  not  Indlcpted  explicitly.) 

When  (2.3)  and  (2,5)  and  higher  functional  derivatives  are  substituted 
Into  (2.4),  we  obtain  a  formally  exact  Infinite  series  for^.  If  tpt  Is 
chosen  to  bo  I,  I.e.  uo  *  0  which  corresponds  to  an  Ideal  gas,  this  series 
becomes  the  usual  virlal  expansion.  In  this  case  Ap(r)  Is  the  Mayer  f 
function,  -1,  and  the  first  functional  derivative  term  becomes  the 


familiar  second  virlal  coefficient  contribution,  since  y«  ■  1  for  an  Ideal 


gAS.  The  second  functionol  derivative  term  vanishes  because 
arc  Jiero  for  on  Idcnl  gas.  Higher  functionol  derivatives  give  the  remain¬ 
ing  terms  of  tlio  virlal  series. 

To  discuss  soft  spheres,  ha/ever,  vra  v/111  use  this  expansion  with  the 
choice  ^0  Pj  Aod  s  -  (p^.  The  resulting  equation  Is  simplified  If 
we  define 

BjCO  -  VdCO^PgCr)  -Oj,(r)J. 

This  function  Is  shavn  In  Figure  Id.  As  can  bo  seen,  It  Is  natural  to  call 
this  the  "blip"  functloa  Then  Fq.  (2.4)  can  be  expressed  as 

( P*/2)Jd£B^{ r) 

+  (p  V2V)rd£>B^(  r , ,)  Bj,(  r ,  s)  (£») 

4  (pVBV)rdj^*B^(r,,)Bj(rs4)Jj^*)(r*)  (2.6) 

4  ... 

The  subsequent  terms  contain  three  or  more  factors  of  B^(r)  (or  p^(r)  '‘p^(r)) 
which  Is  non-zero  only  for  a  small  range  of  values  of  r  (namely, 

for  |r-d|  <  Fd).  Hence  v«  might  expect  that  tlie  contribution  of  the  nth 
functional  derivative  term  Is  of  order  We  have  not  yet  chosen  d  however. 
It  Is  reasonable  to  choose  d  such  that 

yd(»’)r<P,(»’)  -  VdCOl  -  Jdr  Bj(r)  ^  0.  (2.7) 

This  choice  causes  the  first  functional  derivative  term  (which  Is  apparently 
of  order  ^  and  hence  Is  potentially  the  largest  term)  to  vanish  Identically. 
This  Is  Indeed  a  felicitous  choice  for  d,  since  It  makes  the  second  derivative 
terms,  which  nominally  appear  to  bo  of  order  actually  become  of  order 

(see  Appendix  B).  Indeed  It  can  be  shown  that  all  the  subsequent 
terms  are  of  order  5*  or  higher  order.  Thus  we  have  the  result  that 


(2.8) 
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C2(  PiPl'vT's)  ^  ^(p.PJtPd)  (•  + 

when  d  is  chosen  to  satisfy  Eq.  (2.7).  Equation  (2.  C)  provides  a  direct 
connection  between  the  thermodynamic  properties  of  hard  spheres  and  soft 
spheres.  The  associated  hard  sphere  diameter  is  density  and  temperature 
dependent,  and  (as  will  be  seen  below)  has  the  physically  reasonable  property 
that  d  decreases  as  the  temperature  or  the  density  is  increased. 

The  same''blip  function"  expansion  con  be  used  to  obtain  an  expression 
for  y^{r)  using  Eqs.  (2.6)  and  (2.3). 

+  (p'/JV);()r‘6(r,,-r)Bj(rj.)J^(‘)(r‘)  ...)  . 

In  generai  it  is  not  possible  to  choose  d  such  that  the  first  functional 
derivative  correction  term  for  y  vanishes  Identically  for  all  r.  This 
first  correction  appears  to  be  of  order  However,  if  we  use  the  same 
value  of  d  as  v;o  chose  abovu  for  the  free  energy  series.  It  con  be  shown 
(seo  Appendix  B)  that  this  term  is  actually  of  order  for  all  r.  Indeed 
all  subsequent  terms  are  of  order  or  smaller.  Thus  we  have 

9,(r)  ,  yj{r)  ()  ♦  0(5>))  (2.9') 

This  equation  gives  a  relationship  between  the  structure  of  a  soft  sphere 
fluid  and  that  of  a  hard  sphere  fluid. 

Equations  (2.6)  and  (2.9'),  together  with  (2.7),  provide  convenient 
approximate  expressions  for  the  excess  Helmholtz  free  energy  and  pair 
correlation  function  of  a  fluid  with  a  repulsive  Intermolecular  potential. 

To  opply  these  formulas  It  is  necessary  to  have  information  about  the  excess 
Helmholtz  free  energy  and  y(r)  for  the  hard  sphere  fluid.  Fortunately, 


Vcrlot  and  Wcls*^  have  presented  analytic  uy.|»resslons  for  tliesc  functions 
v/hich  accurately  sunxiiarize  the  results  of  ccviiputer  calculations  on  hard- 
sphere  systems. 

In  previous  work  by  the  present  authors  concerning  the  effect  of 
repulsive  forces  on  the  structure  of  the  Lennard-Joncs  fiuid,  Eqs, 

(2. P)  and  (2.9),  v/i thout  the  correction  terms,  and  (2.7)  were  postulated 
(on  the  basis  of  some  intuitive  arguments)  to  be  the  relationships  between 
hard  sphere  and  soft  sphere  fluids.  It  is  now  seen  that  these  postulates 
are  really  the  leading  terms  in  systematic  expansions  in  powers  of  the 
softness  of  the  potent iai.  It  can  also  be  seen  that  the  thermodynamic  re¬ 
lationship  (2.  C)  is  inherently  more  accurate  than  the  structural  one  (2.9). 
One  might  expect  the  former  to  be  accurate  even  for  quite  soft  potentials 
since  the  first  correction  is  of  fourth  order  in  the  softness. 

When  the  excess  energy  and  excess  pressure  arc  obtained  by  differen¬ 
tiating  the  lo-west  order  result  for  the  free  energy,  Eq.  (2.8),  the  errors 
are  of  order  If  these  excess  quantities  arc  calculated  from  the 
radial  distribution  function  using  the  energy  and  virial  equations, 

PAE/N  «  (0p/2)fd£  g(r)  u(r)  ,  (2.10) 

pAp/p  =  -(po/6)  f  d£  g(r)  r  (du(r)/dr)  ,  (2.11) 

the  errors  in  each  will  be  of  the  same  order  as  the  errors  in  g,  namely 
if  the  lowest  order  g  is  used  and  if  the  next  approximation  is  used. 

It  can  be  shown  however  that  the  ratio  of  the  excess  energy  to  the  excess 
pressure  is  the  same  for  both  approximations  for  g  and  Is  the  same  as  the 
ratio  obtained  v/hen  the  energy  and  pressure  are  calculated  from  the  lowest 

approximation  to  the  free  energy.  Since  this  latter  ratio  has  error 
only  of  order  it  foil  ays  that  the  ratio  of  excess  energy  to  excess 
pressure  calculated  from  either  of  the  two  lowest  order  approximations 
for  g  Is  more  accurate  than  the  individual  values. 
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The  blip  function  expansion  v^as  developeci  specifically  to  apply  to 
dense  fluids,  Ha/cver,  it  is  easily  shov/n  that  the  theory  is  exact  In 
the  limit  of  lav  dens i ty.  The  correct  second  virial  coefficient  is 
obtained  from  Eq,  (2.0),  and  the  correct  g(  r)  is  obtained  from  Eq, 
(2.9').  As  a  necessary  corollary  of  this  it  follows  that  the  numerical 
value  of  a  correction  of  given  order  In  F,  depends  on  tiie  density  and 
may  be  significant  at  high  density  although  negligible  at  low  density. 
In  the  next  section  we  examine  the  numerical  predictions  of  the 
theory. 


in.  Nuinffrlcal  Results 


A,  Inverse  Twelfth  Po^/er  Potential 

Kecent  Monte  Carlo  calculations  for  a  "soft  sphere"  fluid  with  an 
inverse  twelfth  power  pair  interaction  potential 

u(r)  -  c(o/r)'*  (3.1) 

provide  accurate  thermodynamic  data  v/ith  v/hich  to  compare  our  theory. 

This  potential  is  sufficiently  "soft"  that  it  provides  a  severe  test 
of  the  accuracy  of  our  metliod. 

There  is  no  temperature- independent  characteristic  length  in  an 
inverse  power  potential.  As  a  result,  the  excess  thermodynamic  func¬ 
tions  have  simple  scaling  properties.  For  example,  for  the  potential 
in  Eq.  (3.1)  it  can  be  shown  that  ^AA/N  (  =  -^/p  )  is  a  function  only 
of  the  single  variable  (bc)'^‘  po’.  This  scaling  law  has  two  immediate 
consequences.  From  the  thermodynamic  equations 

pfip/p  “  -  p  (Ij  (^/p))  (*•*) 

and 

WE/H  =-  5  (1^  (^2/p))  (3.3) 

it  follows  that  PAp/p  and  PAE/N  are  also  functions  of  (pc)'^*  po’  only 
and  are  related  by 

PAp/p  =  4  PaE/N  .  (3.4) 

An  analogous  scaling  law  holds  for  the  radial  distribution  function 
9('’i3iP)  y(>‘:P,p).  These  functions  actually  depend  on  only  two 
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1  iKlcpcivk'iil  Viiriiiblcs,  o  scaIccI  length  which  is  r/o(Pe)'^'*  end  the 
single  thurriodynninlc  vorlablo  (gc)'/*  pa*;  I.e. 

?i(r;lkp)  =  G  (r/oOc)'/'»;  Oc)'/*  po*)  (3.5) 

v/htTC  G  is  u  universal  function  for  the  r”'*  potential. 

As  a  result  of  these  scalings  laws,  Monlc  Carlo  calculations  for 
this  poteiuk.l  need  be  performed  for  only  one  temperature.  The  thermo¬ 
dynamic  prepurlies  and  the  radial  distribution  functions  along  one  isotherm 
con  be  used  to  obtain  the  results  for  any  other  temperature  and  density. 

The  scaling  laws  hold  for  the  exact  thermodynamic  properties  and 
pair  correlation  function.  The  results  of  an  approttimate  theory  might 
not  obey  these  relationships.  The  Rowlinson  theory,  for  example,  does 
not  obey  the  scaling  laws.  For  the  blip  function  expansion,  however, 
it  can  be  slK'wn  that  each  of  the  terms  in  the  series  for  3AA/N  satisfies 
the  scaling  law  exactly,  i.e,  is  a  function  only  of  po*,  pro¬ 

vided  that  the  hard  sphere  diameter  d  is  chosen  according  to  Eq.  (2.7). 
Furthermore,  each  term  in  the  blip  function  expansion  for  y  satisfies 
the  scaling  law  implied  by  Eq.  (3.5). 

We  now  compare  the  results  of  the  theory  with  Monte  Carlo  cal¬ 
culations  and  other  theoretical  methods.  We  will  use  only  the  leading 
term  in  the  blip  function  expansion  for  both  the  free  energy,  Eq.  (2.6), 
and  y(r),  Eq.  (2.9).  The  computational  procedure  is  very 

simple.  Along  one  isotherm,  Eq.  (2.7)  is  solved  at  a  number  of 
densities  to  give  the  associated  hard  sphere  diameter  d(p,p).  (See 
Appendix  C  for  a  discussion  of  the  hard  sphere  information  that  is 
needed  for  this  calculation.)  At  each  density  the  excess  free  energy 
of  the  associated  hard  sphere  system  is  calculated,  and  by  Eq.  (2.6) 
this  is  equal  to  the  excess  free  energy  of  the  soft  spheres  under 
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conolfli;  lilt  Ion,  The  compress  I  bt  11  ty  factor,  fp/p,  Is  obtolncd  from 
Eq,  (3.2).  Note  that  since  d  Is  a  (decreasing)  function  of  the  density, 
the  pressure  of  the  soft  spheres  Is  not  simply  equal  to  that  of  the 
associated  hard  split; re  system.  Instead  v/c  have 


(5.6) 


The  secoii'!  term  In  square  brackets  lov/ers  the  pressure. 

In  Table  I,  v;o  compare  the  results  of  our  theory  v/lth  the  Monte 
Carlo  colculat  Ions, '  * » '  *  the  theory  of  Ro/vl  1  nson' (as  generalized  by 
Barker  and  llenJerscn' *) ,  and  the  variational  theory  of  Barker  and 
Henderson''^  on  the  Isotherm  £c-1.  These  theories  ore  briefly  de¬ 
scribed  above  In  Section  I. 

It  can  be  seen  that  the  Rowllnson  method  gives  reasonably  accurate 
free  energies  and  pressures  at  lower  densities,  but  It  breaks  down  at 
‘high  density.  This  Is  to  be  expected  from  a  theory  v/hich  attempts 
to  dcscrihi!  soft  spheres  v;lth  a  density- 1 nde pendent  effective  hard 
sphere  diameter.  The  Barker-Hcndcrson  variational  method  allows  the 
hard  sphere  diameter  to  be  density  dependent  and  hence  gives  consistently 
better  results  over  a  range  of  density  (except  for  low  densities  v/here 
apparently  the  variational  method  fails  to  give  the  correct  second 
virlol  coefficient).  The  results  of  the  blip  function  expansion  are 
better  still.  They  ore  exact  In  the  limit  of  low  density.  The  free 
energy  values  (over  the  entire  density  range  Indicated)  end  the  pressures 
(except  for  the  highest  density)  are  all  within  a  percent  of  the  Monte 
Carlo  values. 
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Thus  v/c  con  concludo  that  Eqs.  (2.8)  and  (2.7)  provide  a  very 
accurate  rehitionship  between  the  thermodynamic  properties  of  hard 
spheres  and  soft  spheres  with  Inverse  tv^elfth  power  repulsions. 

Since  Monte  Carlo  pair  correlation  functions  have  not  been  pub¬ 
lished  for  this  potential  we  cannot  test  the  structural  relationship, 

Eq.  (2.9'),  directly.  V/c  can  provide  on  Indirect  test,  hov/ever,  by 
using  (2.9'}  to  calculate  virlal  pressures  according  to  (2.11).  The 
virial  pressure  results  are  not  expected  to  bo  as  accurate  as  the  free 
energy  pressures  In  Table  Ib,  since  the  latter  have  errors  of  order 
end  the  former  have  errors  of  order  It  can  be  seen  from  Table  II 
that  the  virlol  pressures  are  accurate  at  low  densities  but  are  In 
error  by  as  much  as  12^  at  the  highest  density. 

In  summary,  Eq.  (2.8)  for  the  free  energy  of  a  soft  sphere  fluid  Is 
very  accurate  for  the  Inverse  twelfth  power  potential,  but  Eq.  (2.9*) 
for  the  radial  distribution  function  Is  very  accurate  only  at  low 
densities  for  this  soft  a  potential.  For  a  harder  potential  both 
equations  are  more  accurate,  as  v^e  shall  see  In  the  fol laving  portion 
of  this  section. 
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B.  The  Potential  Which  has  the  Same  Repulsive  Forces  as  the  Lennard- Jones 
Potential 


In  an  earlier  discussion  of  the  role  of  repulsive  forces  In  deter¬ 
mining  liquid  structure,'*  we  considered  a  fluid  whose  Intcrmolecular 
potential  was  purely  repulsive  and  had  exactly  the  same  repulsive  inter- 
molecular  force  as  the  Lennard-Jones  potential  W|^(r),  For  this  fluid, 
the  pair  Interaction  is 


=  0 


r<2 


1/6 


r>2 


1/6 


o 


a 


where  -c  is  the  depth  of  W|^(r)  at  Its  minimum  value.  The  potentials 
W|^  and  Uj^i^  have  exactly  the  same  repulsive  forces  but  the  latter  has 
no  attractive  forces.  For  the  purpose  of  this  paper,  U|^|^  is  of  interest 
because  it  is  purely  repulsive  and  because  Monte  Carlo  calculations 
of  its  structure  and  thorrmodynamic  properties  have  been  performed.  *»  ” 

The  Monte  Carlo  values  of  the  pressure  and  excess  internal  energy 
are  shown  in  Table  III,  where  they  are  compared  with  results  obtained 
from  Eq.  (2.8),  together  with  (3.3)  and  (3.4).  This  provides  a  test 
of  the  thermodynamic  relationship  between  soft  and  hard  spheres. 

The  pressures  predicted  by  the  generalized  Rowlinson  method  are  also 
shown.  At  high  densities  the  latter  differ  from  the  Monte  Carlo  values 
by  about  6)(. 

The  blip  function  free  energy  pressures  (Eqs.  (2.8)  and  (3.6)) 
agree  with  the  Monte  Carlo  results  to  within  a  percent  for  all  temperatures 
and  densities  considered.  The  difference  is  of  the  same  size  as  the 
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errors  of  the  Monte  Carlo  calculations  and  of  the  hard  sphere  data  that 
arc  needed  for  the  present  calculations. ' *  The  excess  energies  are  also 
very  accurate.  Thus  we  conclude  that  Eq.  (2.8)  Is  a  very  accurate 
thermodynamic  relationship  between  hard  spheres  and  soft  spheres  with 
this  potential. 

To  test  the  structural  relationship,  Eq.  (2.9'),  the  pressure 
and  energy  can  be  calculated  from  Eqs.  (2.11)  and  (2.10).  The  results 
are  shotrin  in  Table  IV.  The  maximum  difference  between  the  blip  function 
expansion  and  Monte  Carlo  results  is  2)S. 

We  can  provide  a  more  direct  test  of  the  structural  ralationship 
in  the  Lennard- Jones  repulsion  system  because  Barker  and  Hendarson 
have  obtained  Monte  Carlo  values  for  g(r).'*  Tha  comparison  is  shown 
in  figure  2  for  the  first  peak  in  the  pair  correlation  function  for  one 
temperature  and  density.  The  major  noticeable  differences  between 
the  Monte  Carlo  and  blip  function  expansion  results  era  that  the  peak 
heights  differ  by  S$b  and  the  depths  of  the  minimum  near  r«1.50  differ 
by  about  .05.  At  other  parts  of  the  curve,  the  difference  is  at  most 
.03.  (A  part  of  this  difference  can  be  explained  by  the  inherent 
errors  of  the  Monta  Carlo  method.'*)  This  represents  remarkebly  close 
agreement  and  provides  a  very  strong  confirmation  of  the  structural 
relationship  (2.9')  between  soft  and  hard  spheres  for  this  hard  a 
repulsive  potential. 

In  summary,  for  this  potential,  which  is  much  harder  than  the 
Inverse  twelfth  power  repulsion,  Eqs.  (2.8)  end  (2.9'),  together  with 
(2.7),  are  very  accurate  aquetlons  for  the  thermodynamic  and  structural 
properties  of  the  fluid. 
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C.  Remarks 

I.  A  limitation  of  the  method.  This  method  relies  on  knowledge 
of  the  equation  of  state  end  pair  correlation  function  of  the  herd 
sphere  fluid.  Since  herd  spheres  undergo  a  phase  transition  to  e 

sol  Id- like  phase  when  pd’M.93,  this  provides  an  upper  limit  on  the 
densities  at  which  the  method  as  vw  have  described  It  can  be  applied. 

However,  at  higher  dansi ties  can  at  least  obtain  the  thermodynamic 
properties  by  e  variational  procedure  similar  to  that  of  barker  end 

Henderson,"^  except  using  the  blip  function  method  rather  than  tha 
general  laad  Rowllnson  method  to  treat  the  hard  sphere  aspects  of  the 
calculation.  Real  liquids  of  course  have  attractive  as  well  as 
repulsive  forces,  end  the  present  considerations  ere  only  the  start¬ 
ing  point  for  more  complete  descriptions  of  liquids.  At  high  den¬ 
sities,  real  liquids  as  well  as  herd  sphere  liquids  solidify,  end 
so  the  density  limitations  of  the  present  mathod  ere  not  Important. 

II.  Density  end  temperature  dependence  of  the  herd  sphere 
diameter.  The  crucial  link  between  herd  end  soft  sphere  fluids  Is 
provided  by  Eq.  (2.7)  which  tells  hew  to  choose  the  essecleted  hard 
sphere  diameter  for  e  soft  potentlel.  Numerical  solution  of  this 
equation  shows  that  the  diemeter  depends  on  density  end  temperature 

In  e  physically  reasonable  way  (see  Tables  II  end  IV  and  Fig.  2  of  Raf.  12). 

As  the  temperature  Is  Increased  (at  constant  density),  the  diameter  decreases. 
This  reflects  the  physical  fact  that  at  higher  temperatures  the  Inter- 
molecular  configurations  of  high  potentlel  energy  become  more  prob- 
dble,  end  these  configurations  correspond  to  smaller  Intermoleculer 
separations.  Also,  as  the  density  Is  Increased  at  constant  temper¬ 
ature,  the  diameter  decreases.  This  reflects  the  physical  fact  that 
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at  higher  densities  the  pressure  Is  greater  and  the  molecules  are  squeezed 
closer  and  closer  together. 
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IV.  Conclusion 

The  purpose  of  this  paper  Is  to  present  and  Justify  a  certain  way  of 
relating  the  hard  sphere  fluid  model  to  more  realistic  soft  sphere  repulsive 
potentials  In  the  equilibrium  statistical  mechanics  of  fluids.  We  have  shown 
that  the  Helmholtz  free  energy  and  the  function  y(r)  =  g(r)  for  a 

soft  sphere  fluid  are  approximately  the  same  as  the  corresponding  quantities 
for  a  hard  sphere  fluid  at  the  same  density,  provided  that  the  hard  core 
diameter  Is  chosen  properly.  A  prescription  Is  given  for  calculating  this 
diameter  which  has  the  physically  reasonable  property  of  being  a  decreasing 
function  of  temperature  and  dansity.  These  results  are  derived  from  well  de¬ 
fined  expansions  of  the  Helmholtz  free  en-rrgy  and  y(r)  In  powers  of  a 
parameter  which  Is  a  measure  of  the  softness  of  the  repulsive  potential  or, 

In  other  words,  the  extent  to  which  the  potential  differs  from  the  hard 
sphere  model  potential.  Calculations  based  on  this  theory  are  not  computa¬ 
tionally  difficult.  They  require  detailed  but  available  knowledge  of  the 
equilibrium  properties  of  hard  sphere  fluids.  The  theory  becomes  exact  In 
the  limit  of  low  density,  but  the  results  remain  useful  and  accurate  up  to 
quite  high  densities,  depending  on  the  softness  of  the  potential.  The  upper 
limit  to  the  density  at  which  the  calculations  can  be  performed  Is  determined 
by  the  highest  density  at  which  we  have  accurate  hard  sphere  radial  distri¬ 
bution  functions.  This  density  is  approximately  that  at  which  the  hard 
sphere  system  undergoes  Its  phase  transition. 

In  this  paper  we  have  restricted  our  attention  to  intermolecular  poten¬ 
tials  which  are  positive  and  repulsive.  In  real  fluids,  however,  the 
intermolecular  forces  are  attractive  for  some  ranges  of  Intermolecular  separa¬ 
tions.  Nevertheless,  at  high  densities  (and  at  high  temperatures  for  all 
densities)  the  structure  of  a  fluid  Is  dominated  by  the  repulsive  forces," 
and  so  an  accurate  theory  of  repulsive  forces  can  provide  a  foundation  for 
an  equilibrium  theory  of  liquids."** 
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Appendtx  A:  Evoluotlon  of  Functional  Derivatives  of  ^ 

The  functional  derivatives  of^v/lth  respect  to  cp(£)  can  be  evaluated  In 
at  least  two  ways.  One  way  Is  to  express  ^In  terms  of  the  configuration  In¬ 
tegral  0  (Eq.  (2.1))  and  take  the  functional  derivative  for  finite  V.  The 
result  can  bo  expressed  In  terms  of  reduced  correlation  functions  In  the  can¬ 
onical  ensemble  for  a  finite  volume.  Then  the  thermodynamic  limit  can  be  taken. 
For  second  and  higher  derivatives,  this  method  Is  tedious  because  In  the  can¬ 
onical  ensemble  the  Ursel I  cluster  functions  do  not  approach  zero  as  the  posi¬ 
tion  coordinates  arc  separated,  but  rather  there  arc  extra  l/V  components  to 
these  functions. 

A  somev/hat  easier  method  Is  to  begin  with  the  thernx>dynamlc  limit  virlal 
(cluster)  series  for  The  functional  derivatives  of  this  Infinite  series  are 
relatively  easy  to  calculate,  and  the  resulting  Infinite  series  can  be  expressed 
In  terms  of  correlation  functions  with  the  aid  of  diagrammatic  manipulations 
and  the  product  theorem  for  generating  functions  of  classes  of  diagrams.'* 

By  using  this  latter  procedure,  Eq.  (2.5)  Is  straightforwardly  obtained. 

The  function  ’)(£'’)  Is  defined  as 

Syto  ••• 


1 


m=o 


where  Is  the  sum  of  contributions  from  all  distinct  labeled  graphs 

with  3  root  points  at  £iiX,ii£'3i  ®nd  m  field  points  at  £4f*i£,^3  that: 

(a)  at  most  one  Mayer  f  bond  (Be'^^-I),  connects  each  pair  of  points,  (b)  no 
f  bond  connects  points  I  and  2,  (c)  no  f  bond  connects  1  and  3,  (d)  the  dia¬ 
gram  becomes  a  star  (multiply  connected)  If  fu  and  fu  are  added,  (e)  the 
diagram  does  not  become  disconnected  If  points  I  and  2  are  removed;  and  (f) 
the  diagram  does  not  become  disconnected  If  points  1  and  3  are  removed.  It 
can  be  shown  that  this  Is  precisely  the  series  for  the  function  on  the  right 
hand  side  of  the  following  equation 


(A.1) 


j( ’)(£’)  =  h(r2  5)  +  g(r2,) 


v/hero  h(r)  -  g(r)  -I,  and  g^’)  is  the  three  particle  correlation  function. 

This  final  expression  now  contains  no  reference  to  diagrams.  (The  same 
result  can  also  be  obtained  by  the  first  method  described  in  this  Appendix.) 

If  one  makes  the  Kirkwood  superposition  approx imat ion, we  would  have  simply 

'  h(r.O.  (A.*) 

Similarly  the  function  Is 

m=o 

where  the  integrand  is  the  sum  of  contributions  from  ail  distinct  labeled 

graphs  with  4  root  points  (1, 2,3,4)  and  m  field  points  (5,  ...,  iih4)  such  that: 

(a)  at  most  one  f  bond  connects  each  pair  of  points,  (b)  no  f  bond  connects  1  and  2 

(c)  no  f  bond  connects  3  and  4,  (d)  the  diagram  becomes  a  star  if  f|i  ^s« 
are  added,  (e)  the  diagram  does  not  become  disconnected  if  points  1  and  2 

are  removed,  and  (f)  the  diagram  does  not  become  disconnected  if  points  3  and 

4  are  removed.  An  expression  for  j(*}  in  terms  of  four,  three,  and  two 

particle  distribution  functions  can  also  be  obtained.  If,  in  addition,  we  make 

the  superposition  approximation,  we  have 

j(*)(£*)  »  h(r,,)h(rj4)  +  h(ru)h(r23)  +  h(r,,)h(ri»)h(r,4) 

+  h(ris)h(r,4)h(r,4)  +  h(r,,)h(ri4)h(r,,) 

+  h(ri4)h(ris)h(r,4)  +  h(ris)h(ri4)h(rt5)h(rt4).  (a.3) 

This  result  can  also  be  obtained  by  using  the  diagrammatic  expansion 

for  and  summing  only  those  graphs  which  are  consistent  with  the 

to 


superposition  approximation. 
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Appcndix  B:  Estimating  the  ^  Dependence  of  Various  Integrals 

The  effective  hard  core  diameter  d  is  chosen  according  to  Eq.  (2.7). 

As  a  result  of  this  choice,  many  integrals  arc  automatically  of  higher 
order  in  §  than  one  might  ordinarily  expect.  In  this  Appendix  we  discuss  the 
reason  for  this. 

By  converting  the  integral  in  Eq.  (2.7)  into  a  one  dimensional  integral 
we  find 

•» 

J  r'dr  Bj(r)  =  0.  (B.  1) 

o 

Consider  an  integral  of  the  following  type 

f  r*dr  Bj(r)  F(r;s)  (B.2) 

o 

where  F  Is  any  smooth  function  of  r  and  may  depend  on  some  other  variables, 
denoted  by  s,  which  are  held  constant  In  the  r  Integration.  Since  B|j(r) 

Is  effectively  non-iero  for  |r-d|  <  5d  only,  we  might  expect  this  integral 
to  be  of  first  order  In  ?.  However,  because  of  Eq,  (B.1),  wa  have 


J  r'dr  B.(r)  F(r;s)  =  J  r'dr  B.(r)[F(r;s)  -  F(d;s)) 
o  o 

If  F(r;s)  Is  a  slowly  varying  function  of  r  (I.e. ,  If  F  Is  differentiable) 
then  the  quantity  In  square  brackets  Is  of  order  ^  and  the  Integral  Is  of 
order  $',  not 


To  apply  this  rasult  let  us  no#  consider  the  Integral  containing 
In  the  right  hand  side  of  Eq,  (2,6),  Fix  jji  and  and  consider  doing 
the  £t  Integration.  The  above  result  cannot  be  applied  directly  because 
Jd^’^(r’)  Is  not  a  continuous  function  of  |rit|  at  fixed  n,  rs,  and  fixad  direction 
for  ri'  (see  Eqs.  (A,  1)  and  (A, 2)).  However,  we  first  Integrate  over  all 
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the  tingles  ossoclotccl  with  rit,  holding  Iritl  fixed.  The  discontinuities 
arc  smoothed  out  and  nov  the  integration  over  l^itl  gives  a  result  of  order 
F,*.  The  remaining  integral  over  Irisj  gives  another  factor  of  5**  Thus 
the  complete  integral  is  of  order 

This  method  can  be  employed  to  verify  all  the  statements  given  in  the 
text  concerning  the  §  dependence  of  various  integrals. 
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Append  lx  C:  Information  about  the  Hard  Sphere  Fluid 

To  perform  calculations  using  the  present  method,  It  Is  necessary 
to  kna-j  the  equation  of  state  for  hard  sphere  fluids  as  well  as  the 
function  yj( r) . 

The  equation  cT  state  which  we  use  Is  the  one  suggested  by 
Carnahan  and  Starling, 

3Pj/P  (HiH  V- V)/(1-Ti)*  , 

where  q  =  TTpd^/6  and  d  Is  the  hard  sphere  diameter.  This  formula 
summarizes  the  available  molecular  dynamics  and  Monte  Carlo  results 
within  the  statistical  accuracy  of  these  computer  calculations.  By 
Integrating  with  respect  to  density,  the  excess  Helmholtz  free  energy 
corresponding  to  this  equation  of  state  can  be  obtained: 

“  -ti(A-3ti)/(1  -Ti)*  . 

The  y^(r)  function  Is  also  needed  for  all  values  of  r>d  as  well 
as  for  r  slightly  less  than  d,  I.e.  Just  Inside  the  hard  core.  For 
r>d,  y^(r)  «  9^(r),  which  can  be  obtained  from  Monte  Carlo  and  mole¬ 
cular  dynamics  calculations.  We  use  the  empirical  formula  of  Varlet 
and  Weis*,  which  expresses  as  the  sum  of  two  parts: 

*  9i(r)  +gt(r),  r>d  ,  (C.l) 

where  gi(r)  Is  simply  related  to  the  Werthelm-Thel le*  analytic  solu¬ 
tion  of  the  hard-sphere  Percus-Yevick  equetlon'  and  gt(r)  Is  a  small 
correction  (also  In  convenient  analytic  form).  Verlet  and  Weis  state 
that  g^(r)  obtained  In  this  simple  way  differs  from  their  Honte  Carlo 
results  by  at  most  0.03  and  estimate  the  statistical  erior  to  be  approximately 
0,01.  To  obtain  y^  for  r  near  d  Inside  the  core,  we  assume  we  can 
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uso  the  same  formula.  Inside  the  core,  the  Porcus-Yevlck  1$  known 
analytically,  and  the  value  of  0}  Is  obtained  by  extrapolation.  We 
adopted  the  procedure  of  extrapolating  the  logarithm  of  gj  quedretical ly 
Into  the  core.  This  procedure  is  clearly  arbitrary;  however,  comperlson 
of  different  extrapolation  methods  indicates  that  an  error  of  et  most 
a  few  parts  in  10*  are  introduced  into  the  calculated  d  by  the  arbi¬ 
trariness  of  tlie  extrapolation.  For  exemple,  the  values  for  d  given 
in  Table  IV  differ  from  those  calculated  by  Verlet  end  Weis  by  e  differ¬ 
ent  extrapolation  procedure  by  at  most  .0002.  The  error  Is  so  small 
because  d  depends  on  y^  inside  the  core  only  for  values  of  r  which  are 
very  close  tod  end  any  reasonably  extrepolet Ion  procedure  gives  much 
the  same  results.  This  error  makes  e  negligible  difference  for  the 
thermodynamic  properties  celculeted  from  the  free  energy.  The  extra¬ 
polation  procedure  can  obviously  hove  e  small  effect  on  the  value  of 
g(r)  for  r<d  when  It  Is  calculated  from  Eq.  (2.9*),  especially  for 
very  soft  potentials.  This  effect  is  minimized,  however,  by  the  ex¬ 
ponential  factor  which  approaches  zero  rapidly  es  r  Is  decreased  from 
d.  Nevertheless,  this  small  effect  on  g( r)  could  heve  en  effect  on 
the  virlal  pressure  and  energy  calculeted  from  Eqs.  (2.11)  and  (2.10) 
because  of  the  sensitivity  of  these  equations  to  slight  errors  In  g. 

We  estimate  that  for  tha  Inverse  twelfth  power  potential,  the  virlel 
pressures  might  contain  an  error  of  the  order  of  et  most  S)(  et  tha 
highest  density  (and  much  less  at  lower  densities),  while  for  the  much 
harder  Lennard-Jones  repulsive  potential,  tha  error  In  the  virlel 
pressure  due  to  the  uncertainty  In  the  extrapoletlon  Is  negligible 
compared  with  errors  erlsing  from  Inaccuracies  In  our  g|||(r)  for  r  outside 
the  core.  These  latter  errors  are  of  tha  order  of  a  percent.'* 
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Table  lo.  Reduced  Excess  Free  Energy  per  Particle  for  a  Fluid  of 

Molecules  Interacting  Through  an  Inverse  Twelfth  Power  Repulsion*  — 
A  Test  of  the  Thermodynamic  Relationship  between  lie rd-Sphe remand 
Soft-Spheres . 


/p  =  3AA/N 


po* 

Monte  Carlo^ 

Monte  Carlo^ 

This  Work** 

8arkor-llcnderson 

Variational 

Ganerallzed 
Rowl 1 nson 

.UU 

.40 

.40 

.40 

.39 

.39 

.2828  , 

.91 

.91 

.91 

.89 

.89 

.4243 

1.53 

1.53 

1.54 

1.54 

1.54 

.5657 

2.32 

2.33 

2.33 

2.37 

2.42 

.7071 

3.33 

3.34 

3.34 

3.43 

3.68 

.8485 

4.60 

4.61 

4.65 

6.77 

5.58 

a.  These  celculallons  are  for  the  Isotherm  3c  1. 

b.  Reference  U. 

c.  Reference  15. 

.  Calculated  from  Eqs.  (2.8)  and  (2.7). 
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Tablc  !b.  Pressure  of  a  Fluid  of  Molecules  Interacting  Through  an 
Inverse  Twelfth  Power  Repulsion^  --  A  Test  of  the  Therfflodynamlc 
Relationship  ilctviecn  Hard-Spheres  and  Soft-Spheres. 

Pp/p 


O 

a. 

Monte  Carlo^ 

Monte  Carlo® 

This  Work** 

Barker-Hendersen 

Variational 

General  I  zed 
Rowl  in'son 

.1414 

1.45 

1.45 

1.45 

1.36 

1.44 

.2828 

2. 12 

2. 1? 

2.12 

2.11 

2.13 

.4243 

3.  to 

3.12 

3.12 

3.18 

3.25 

.5657 

4.56 

4.58 

4.57 

4.72 

5.16 

.7071 

6.64 

6.66 

6.71 

6.93 

8.57 

.8485 

9,46 

9.56 

9.89 

9.98 

15.18 

a.  These  calculations  are  for  the  Isotherm  &  1. 

b.  Reference  14. 

c.  Reference  15. 

d.  Calculated  from  Eqs.  (2.8),  (2.7)  and  (3.2). 


Toblo  II.  Viriol  Pressure  end  Assoctcitccl  Herd  Sphere  Diameter  of  a 
Fluid  of  Molecules  Interacting  Through  an  Inverse  Twelfth  Pov/er 
Repulsion'^  —  A  Test  of  the  Structural  Relationship  Between  Hard- 
Spheres  and  Soft-Spheres. 


^p/p 


po’ 

Monte  Carlo*^ 

This  Work^ 

d/o** 

.1414 

1.45 

1.45 

1 . 0663 

.2828 

2.12 

2.16 

1.0614 

.4243 

3.11 

3.23 

1.0552 

.5657 

4.57 

4.85 

1.0477 

.7071 

6.65 

7.23 

1.0390 

.8485 

9.51 

10.60 

1.0295 

a.  These  calculations  are  for  the  isotherm  ^6  =  1. 

b.  The  average  of  the  two  Monte  Carlo  values  (see  Table  Ib). 

c.  Calculated  from  Eqs.  (2.9*),  (2.7),  and  (2.11). 

d.  The  associated  hard  sphere  diameter,  in  units  of  a,  calculated 
from  Eq.  (2.7). 
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Tablo  HI.  Pressure  and  Reduced  Excess  Internal  Energy  per  Particle 
for  a  Fluid  of  Molecules  Interacting  Through  the  Repulsive  Part  of 
the  Lennard-Joncs  Potential  —  A  Test  of  the  Thermodynamic  Relation¬ 
ship  bctv;ccn  Hard-Spheres  and  Soft-Spheres, 


3p/p  3AE/N 


o 

Q. 

T^® 

Honte  Carlo^ 

This  Work*^ 

General i zed 
Rowl i nson 

Honte  Carlo^ 

This  Work** 

.85 

2.81 

6.92 

6.92 

7.38 

.84 

.84 

.84 

.75 

10.23 

10.33 

10.99 

.97 

.98 

.65 

1.35 

4.89 

4.94 

5.06 

.46 

.46 

.40 

1.35 

2.53 

2.54 

2.55 

.17 

.18 

.85 

.72 

10.87 

10.83 

11.55 

1.03 

.50 

1.36 

3.28 

3.29 

3.31 

.27 

a.  The  reduced  temperature,  kT/e. 

b.  The  first  four  states  are  from  Ref.  6  and  the  last  two  are  from 
Ref. 17. 

c.  Calculated  from  Eqs.  (2.8),  (2.7),  and  (3.2). 

d.  Calculated  from  Eqs.  (2.8),  (2.7),  and  (3.3). 
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Tablc  IV.  V trial  Pressure,  Reduced  Excess  Internal  Energy  Per  Particle, 
end  Associated  Hard  Sphere  Diameter  for  a  Fluid  of  Molecules  Interacting 
Through  the  Repulsive  Part  of  the  Lennard-Jones  Potential  —  A  Test 
of  the  Structural  Relationship  between  Soft-Spheres  and  Hard-Spheres. 


^p/p  3AE/N 


po* 

T^a 

Monte  Carlo^ 

This  Work® 

Monte  Carlo^ 

This  Work** 

d/o* 

.85 

2.81 

6. 92 

7.04 

.84 

.86 

.9699 

.84 

.75 

10.23 

10.40 

.97 

.99 

1.0224 

.65 

1.35 

4.89 

4.99 

.46 

.47 

1.0029 

.40 

1.35 

2.53 

2.55 

.17 

•  V® 

1  1.0048 

.85 

.72 

10.87 

10. 89 

1.04 

1.0237 

.50 

1.36 

3.28 

3.31 

.27 

1.0039 

a,b.  See  Footnotes  of  Table  III. 

e.  Calculated  from  Ef s.  (2.9*),  (2.7),  and  (2.11). 

d.  Calculated  from  Eqs.  (2.9*),  (2.7),  and  (2.10). 

a.  The  associated  hard  sphere  diameter,  in  units  of  a,  calculated  from 
Ed.  (E.7). 
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Figure  1. 


Figure  2. 


Figure  Captions 


Schematic  plots  of  sonK2  functions  considered  in  the  blip 
oxpans ion: 

a.  The  exponential  of  3  times  the  negative  of  a  soft 
sphere  potential. 

b.  The  exponential  of  3  times  the  negative  of  a  hard 
sphere  potential. 

c.  The  difference  between  (a)  and  (b)  showing  the 
significance  of  the  softness  parameter 

d.  The  biip  function,  Bj(r).  According  to  Eq.  (2.7), 

the  diameter  d  is  chosen  to  make  the  net  areas 
under  r^  B^(r)  equal  to  zero. 

The  radial  distribution  function  for  a  dense  fluid  of 
molecules  interacting  through  the  repulsive  part  of  the 
Lennard-Jonos  potential  —  A  test  of  the  structural  relation¬ 
ship  between  hard-spheres  and  soft-spheres.  The  density 
is  po’=.8S.  The  temperature  is  kT/eB.72,  The  dots  are 
the  Monte  Carlo  results  of  Barker  and  Henderson,  Ref.  17. 

The  smooth  curve  was  calculated  from  Eqs.  (2.9')  and  (2.7). 


